Spin-independent effective mass in a valley-degenerate electron system 
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In a generic spin-polarized Fermi liquid, the masses of spin-up and spin-down electrons are ex- 
pected to be different and to depend on the degree of polarization. This expectation is not confirmed 
by the experiments on two-dimensional heterostructures. We consider a model of an ./V-fold degen- 
erate electron gas. It is shown that in the large-JV limit, the mass is enhanced via a polaronic 
mechanism of emission/absorption of virtual plasmons. As plasmons are classical collective excita- 
tions, the resulting mass does not depend on N, and thus on polarization, to the leading order in 
1/N. We evaluate the 1/N corrections and show that they are small even for N — 2. 

PACS numbers: 73.21.-b, 71. 10. Ay, 71.18.+y, 71.10.Ca 
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The observation of an apparent metal-insulator transi- 
tion in high-mobility Si metal-oxide-semiconductor-field- 
effect-transistors (MOSFET's) challenged the scaling 
theory of localization 0], which predicts that a two- 
dimensional (2D) system undergoes only a continuous 
crossover between weak and strong localization regimes. 
Although there has been a substantial progress in under- 
standing of transport and thermodynamic properties of 
MOSFET's and other heterostructures 00, the origin 
of the observed phenomena is still a subject of discussion. 
Although a conventional (dirty) Fermi-liquid (FL) theory 
00 can account for many observed effects at least qual- 
itatively and, in some cases, quantitatively, there is also 
a number of non-FL scenarios for the anomalous metallic 
state [ZllSl- On the experimental side, the main argument 
for the FL-nature of the metallic state is the observation 
of quite conventional Shubnikov-de Haas (ShdH) oscil- 
lations 0, 0] , which implies an existence of well-defined 
quasiparticles albeit with the renormalized effective mass 
to* and spin susceptibility xt- The ShdH and magne- 
toresistance experiments show that at low densities both 
77i* and xt are significantly enhanced compared to their 
band values and, according to some studies 0, [ljj , 
even diverge at the resistive transition point. 

Although none drastically non-FL features of the 
metallic state have been found in ShdH measurements 
as of now, there is one very intriguing observation which 
does seem to present a challenge for the FL theory, at 
least in its conventional formulation. Namely, in all stud- 
ies when the spin and orbital degrees of freedom were 
controlled independently by applying a tilted magnetic 
field, the effective masses, m* and m|, and Dingle tem- 
peratures (impurity scattering rates), Trjj and Tjji, of 
spin-up and -down electrons, were found to be almost 
the same. Moreover, to* in MOSFETs 0,0] was found 
to be independent of the spin polarization, whereas To 
was shown to depend on the polarization only weakly. In 
n-GaAs, the effective mass was found to depend on the 
parallel magnetic field however, this behavior was 
attributed to the coupling between the in- and out-of- 
plane degrees of freedom (Stern effect [13), which is to 



be expected in systems with wider quantum wells. Given 
that the Stern effect is subtracted off, the resulting de- 
pendence of to* on the polarization is likely to be weak. 

Why is this strange? Polarization is expected to lead 
to two effects: the spin-splitting of the effective mass, i.e., 
77i| =/= m* , and dependences of both mi and m* on the 
polarization. The first effect can be understood by con- 
sidering a partially spin-polarized FL as a two-component 
system. As the densities of the components are different, 
the corresponding couplings describing the interactions 
between the same and opposite spins are also different; 
hence a priori the mass renormalizations should also be 
different. That the masses should depend on polarization 
can be seen from considering two limiting cases: of zero- 
and full polarization. At fixed density n, the Fermi en- 
ergy is doubled by fully polarizing the 2D system, hence 
the ratio of the Coulomb to Fermi energy g = e 2 \fim/ Ep 
differs by a factor of 2 between the cases of zero and full 
polarization. The experiment shows that the mass does 
depend on the density; however, if g is the only dimen- 
sionless parameter that determines the mass renormaliza- 
tion, the same effect can be achieved by either varying n 
or by varying Ep via polarization at fixed n. Also, dif- 
ferent Fermi velocities should result in different impurity 
scattering times for spin-up and -down electrons; hence 
the Dingle temperatures are also expected to be different. 
However, this is not what the experiment shows. 

The qualitative arguments given above can be verified 
in a number of ways. Back in 1971, Overhauser pre- 
dicted the spin-splitting and polarization dependence of 
to* within the RPA approximation for the 3D case [lflj . 
Repeating the calculation in 2D gives a similar result: 

m* n /m = 1 + (r s /%/27r) lnr s T (r s £/2\/27r) lnr s , (1) 

(n-f — n{)/ (rij +njj <C 1 is the polariza- 



In the fully-polarized regime 



where £ 

tion and r s = me' 2 / 
(£ = 1) , the spin-down electrons disappear, whereas the 
renormalization of mt is by a factor of \/2 smaller than 
for £ = 0. This argument can be generalized for a (par- 
tially) spin-polarized FL [J^, where the Landau inter- 
action function has three independent components: 
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f^, and /Tl — ylT _ xhe Galilean invariance then gives 
m/mX = 1 - Fp - (k Fl /k F ^)Fl l ; 



m/m* = 1-Fp - {k F] /k Fi )F : 

where F[° = mj d6>cos(9/ iJ (0) /(2vr) 2 , with i,j =f,L 
Again, in general, m* 7^ toT. In addition, the spin- 
splitting and polarization dependence of to* are also ob- 
tained within the Gutzwiller approximation for the Hub- 
bard model 01 ( m this case, mass-splitting disappears 
at half- filling but the polarization dependence survives). 

Absence of the polarization dependence of the effec- 
tive mass suggests that m* is renormalized via the in- 
teraction with some classical degree of freedom, which 
is not affected by the quantum degeneracy of the elec- 
tron states. In this paper, we show such a mechanism 
may be provided by the interaction with (virtual) plas- 
mons which dominate the mass renormalization beyond 
the weak-coupling regime. To this end, we turn to a 
model of a Coulomb gas with large degeneracy A, con- 
sidered previously in Refs. [l^, EH- This model is rele- 
vant, first of all, to valley-degenerate systems, such as the 
(001) surface of a Si MOSFET, where N = 4 (two val- 
leys and two spin projections). As the valley degeneracy 
plays a very important role in the dirty FL theory yj |g 
it is important to elucidate its role for the properties of 
a clean FL. However, the 1/N expansion turns out to be 
converging reasonably fast even for a non- valley degener- 
ate system (N = 2) and, as such, it provides a simple yet 
non-trivial way of going beyond the weak-coupling limit 
for not too strong Coulomb interaction. 

For a 2D A-fold degenerate Coulomb gas, the Fermi 
momentum is scaled down by a factor of TV -1 / 2 (since one 
has to distribute the same number of electrons among A 
isospin flavors), whereas the inverse screening radius (k), 
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FIG. 1: Excitation spectrum for an iV-fold degenerate 2D 
Coulomb gas in the strong-screening regime (r s N 3/2 > 1). 
The plasmon dispersion crosses over from the ^fq to q 2 form at 

q ~ go ~ rl^n 1 / 2 3> k F . Processes with momentum and en- 
ergy transfers in the shaded oval (q ~ qo and lo ~ q 2 , /m) dom- 
inate the mass enhancement. The plasmon spectrum merges 
with the continuum at q = qi ~ rl^ 2 N 1 ^ 4 n 1 ^ 2 ^> qo- 



proportional to the density of states, is scaled up by a 
factor of A. The ratio a = n/k F = r s N 3 ^ 2 /2 controls 
the crossover between the regimes of weak (a< 1) and 
strong (a 3> 1) screening. For N 3> 1, both of these 
regimes are compatible with the condition r s <C 1 which 
guarantees that the screening cloud includes many elec- 
trons, so that the mean-field theory is applicable. For 
a«l, the screening radius nT 1 = cT 1 ^ 1 is larger than 
the Fermi wavelength. [This case also includes the usual 
RPA scheme for A — 2-see Eq. Q.] The mass renor- 
malization is mostly due to elastic scattering within the 
particle-hole continuum with momentum transfers q ~ n, 
whereas the interaction with plasmons is small. In this 
regime, the mass depends on total degeneracy (A) and is 
thus strongly affected by polarization. Also, as scatter- 
ing is mostly by small angles, m* < to. For a»l, the 
effective screening radius q^ 1 = (2a)~ 1 ^ 3 kp 1 is smaller 
than the Fermi wavelength (but still larger than the dis- 
tance between electrons); hence, scattering is isotropic 
(s-wave). The particle-hole continuum contribution to 
to* is greatly reduced for s-wave scattering, whereas the 
interaction with virtual plasmons now plays a dominant 
role. As the plasmon is a classical collective mode, it is 
not affected by a change in A. Consequently, the lead- 
ing term in the A -1 expansion for to* does not depend 
on A, whereas the next-to-leading term happens to be 
numerically small. 

The effective mass is found from the self-energy via the 
usual relation (valid for a small renormalization) 



to* /to = 1 



/ 0AS fc (g) | dA£ fc ( £ ) 
V de k d(ie) 



k— >k F ,e— >0 



where AEfe(e) = Sfe(e) — £fc F (0). It is convenient to 
separate AEfe(e) into the static and dynamic parts as 

A£ fe ( £ HA£f(e) + A£f», (2) 
where the static part for e k = (k 2 — k F )/2m — > is 

£ W^ Vq{0) [Gk+q {£ + w) - GkF+q (£)] 



AEj?(e) 



2tt 

., \2 £k I d0COS&V 2k FS ine/2 
( Z7r ) Jo 



(0) 



with G k 1 (e) — ie — ek and 



(3) 



(4) 



V q (w) = [<?/27re 2 - II g (w)] 1 . 
The dynamic part is 

^fW - /5^iKM-v.(o)] 

x [Gk+ q (e + u) - G kF+q (e)] . (5) 
In what follows, we will need the following two forms 
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of the polarization bubble 



/rip r rl 2 k 
V / -7-2G k ( £ )G k+q ( e + w ) (6) 
2tt J (2ir) 



(mN/2ir) (l - \oj\/y/u, 2 + v 2 F q 2 ) , for g < k F ; 
2ne q / (e 2 + lj 2 ) , for q > k F , 



where v F — y Aim/m 2 N and e q = q 2 /2m. 

In the weak-screening regime, AS|* (e) [Eq. (0)]m* 
gives the main contribution to m*. To logarithmic ac- 
curacy, TO */m = 1 + (r s VN /2tt) In (r s iV 3 / 2 ) + 0(r s ) in 

this regime. [For N = 2 and £ = 0, this reduces back to 
Eq. (UJ]. In this regime, the plasmon contribution to m* 
is a subleading, C(r s )-term. 

Now we turn to the strong-screening regime. The 
static screened potential in Eq. is evaluated for q = 
2k F su\0/2 < 2k F . In this range, V q (0) = 2ire 2 /{q + n) 
is of the same form as in the weak-screening regime but 
now Vjj(0) depends on q only weakly because q <C K. 
Consequently, the angular averaging in Eq. © renders 
the static contribution to to* small: (to*/to— l) st = 
8/3irNa. Using the large-g form of II in Eq. one 
obtains V q (0) = 2ne 2 q 2 / (q 3 + q%) for q > k F , where 

1/3 

<7o = (2a) fcp > fcj? is the inverse screening radius in 
this regime hj. The main contribution to m* comes 
from the region of large q and w in Eq. JSJ, i.e., from 
the plasmon region. In the strong-screening regime, the 

plasmon dispersion is given by lu p — J e 2 + 2ire 2 nq/m. 



The crossover between the ^fq and q 2 behaviors occurs at 
q r-j qo. The plasmon runs into the continuum at q ~ q\ = 
k F {a/2) 1 / 2 q . Most importantly, being the classical 
collective mode, plasmon is not affected by a change in N. 
The mass renormalization can be estimated as follows. 
Typical momenta and energy transfers are of the order 
of q and e qo , respectively; thus V qo {s qo ) ~ e 2 /go, an d 



G 



e q( \ Combining these estimates together, 



one finds that (to*/to- l) dyn ~ J d 2 q J duV q G 



2/3 

which is larger than the static contribution by a 5 / 3 ^> 1. 
To perform an actual calculation, we notice that the plas- 
mon contribution from the region of large q to the effec- 
tive mass can be written as 



to*/to=1H — / (fegRes 
Jo 



(iu - e q f 



(7) 



where only the poles of V q (uj) were taken into account, 
and where we have used the expansion £k+q = £fc + 
v F q cos 9(1 + ek/2E F ) + e q . Substituting the large-g form 
of II [Eq. ©] into V q (uo) in Eq. Q, one arrives at the 
result of Ref. [19| for the leading 1/N term in to* 



to* /to = 1 + Cr 2 J 3 , 
where C = r(l/3)r(l/6)/60V^ ~ 0. 14. 



(8) 



Corrections to the leading term are obtained by includ- 
ing (a) interaction corrections to the bubble [Fig. Ufa)], 

(b) vertex correction to the self-energy [Fig. Hfb)], and 

(c) corrections to the polarization bubble from the small- 
q region. Estimating the diagrams in Fig. |2Ia,b) in the 
same way as for the leading term, we find that both 
(a) and (b) contribute iV-independent, r^ 3 corrections to 
Eq. |JSJ. We have verified by an explicit calculation that 
these estimates do hold. Next, we consider correction (c) 
and show that it gives the next-to-leading term in the 
1/N expansion. 

The I /q correction to the large-g form of the bubble 
[Eq. ®] is 



4n 2 7r (3c^ 2 - el)e 



1> 9 



mN (u 2 + e 2 ) 3 



(9) 



At the plasmon pole {to 2 = —ui 2 ) and for q ~ qo, the 
relative correction \8H q (u)) /U q (u))\ ~ 1/a 2 / 3 , hence one 
can expect the next-to-leading term in the mass to be 
of order r 2 J Z /a 2 ! 3 ~ 1/N. Indeed, a correction to the 
bubble © shifts the position of the plasmon-pole from 
w 2 to w 2 +A 2 , where A 2 = 87r 2 ne 2 (3r 2 + l) /Nmqr 4 and 

r = \ll + (qo/q) 3 - Substituting this result into Eq. J7|), 
and evaluating the g-integral to log-accuracy (the upper 
limit is determined by q ~ q± , corresponding to the region 
where the plasmon runs into the continuum), we obtain 
to* within the next-to-leading order in 1/N as 



1 



to*/to= l+0.14r 2/3 +— — log (r s N 3/2 )+0 



12N 



1 



r s N 5 / 2 

(10) 

where the last term is the static contribution of the con- 
tinuum. We see that the 1/N expansion generates the 
series in powers of (rsN 3 ^ 2 )^ 1 . 

Now we apply our main result, Eq. (|1(J|) . to real sys- 
tems. [In what follows, we neglect the last term in 
Eq. HI Up .] First of all, due to a small numerical coefficient 
in the leading term in Eq. I|1U|) . the actual constraint on 
r s being small is rather soft: a two-fold enhancement 
of the mass occurs only for r s « 20, hence smaller val- 
ues of r s still allows for a reasonable description within 
the mean-field theory. Eq. (|10l) agrees well with the ob- 
served dependence of m* (r s ) for Si MOSFETs in the 
range r s = 2 — 6; for larger r s , the theoretical value 
of to* falls below the experimental one. In the inter- 
val 2 < r s < 6, the 1/N term in Eq. (dJjJ is not that 
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b) 




FIG. 2: a: corrections to the bubble; b: vertex correction to 
the self-energy. 
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FIG. 3: Change in the effective mass under full spin polar- 
ization [cf. Eq. Illl l. as a function of r 3 . Inset: polarization 
dependence of the effective mass for r s — 2, 3, 4, 5. 



small: it constitutes 18-26 % and 26-32% of the leading 
term for N = 4 and N = 2, correspondingly. However, 
the relative change in m* due to full spin polarization 
(N -> N/2) 



= 2 x 



m*(N/2) - m* (N) 



x 100%, 



(11) 



TOavg m*{N/2) +m*(N) 

is small. Am/m avg as a function of r s is shown in Fig. 
for N — 4 and N = 2. In both cases, these changes are 
less than 3 %, which is likely to be below the experimental 
error in the measured mass. At finite polarization, the 
result in Eq. l|lUfl changes to 



-ii = l + 0.14r? /3 + 



12N 



log 



,N 3 / 



2 1 



(12) 



Notice that although an explicit polarization dependence 
does occur in the second term, there is no spin-splitting 
of the masses to this order in X/N. Eq. I|12|) is valid as 
long as there are still many spin-down electrons within 
the screening radius or, equivalently, 1 — £ 3> r 2 J z ~ 
(m*/m— 1). Fig. |21 shows that the effective mass remains 
essentially constant in the whole range of £, which is in 
agreement with the experiment 12] . 

To leading order in l/N, the renormalization of \* 
is entirely due to that in m*, so that g* = x*/ m * re- 
mains unrenormalized 0|. We found that this remains 
true up to the next-to-leading term in X/N. This result 
is in qualitative agreement with the experiments on Si 
MOSFETs. However, recent experiment on AlAs system 
shows that the g* factor is affected by lifting the valley 
degeneracy More work is required to attribute this 
behavior to a many-body effect. 

Now, we comment briefly on the impurity scattering 
rate in the large- N limit. In the strong-screening regime, 
the screening radius (cfo -1 ) is much shorter than the Fermi 
wavelength. Therefore, scattering even on charged impu- 
rities is in the s-wave regime. We assume that the main 



role is played by impurities within the 2D layer. Due to a 
peculiarity of 2D scattering 21], the scale of the scatter- 
ing cross-section section is set by the wavelength (rather 
than by the impurity size a ~ Qo 1 ) and depends on a 
weakly: A ~ fc^/ln 2 (kpa). Consequently, the scatter- 
ing rate 1/r = riiVpA, where rii is the concentration of 
impurities, has only a weak dependence on the polar- 
ization (via Uf under the logarithm). Thus 1/r (Dingle 
temperature) for spin-up and down-electrons are close 
to each other. Notice that both ShdH and weak-field 
Hall effect [2^ show that 1/r, while being the same for 
spin-up and spin-down electrons, increases strongly with 
r s . Within our model, this can only be explained by an 
increase in the number of scatterers rij with decreasing 
electron density-not an improbable scenario for Si MOS- 
FETs. 

Finally, we observe that as the mass is renormalized by 
plasmons with large q, the behavior of the plasmon spec- 
trum at small q (gapped or gapless) is irrelevant. Conse- 
quently, in 3D the mass is renormalized in a similar way: 
m*/m = 1 + C3Dfs^ 4 , where C^d does not depend on N. 
Therefore, a finite thickness of the 2D layer should not 
affect the (approximate) spin-independence of the mass. 
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